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1. Introduction 

A large number of speculative theoretical ideas suggest the existence of new scalar fields, as 
reviewed in [1—3]. A sufficiently light scalar will produce a new macroscopic force. Such forces 
may lead to violations of the equivalence principle (EP) and the inverse square law (ISL). For 
recent discussions of experimental searches for violation of the EP see ref. [4-8], while recent 
reviews of searches for ISL violation may be found in [3,4,9]. 

Recently, Khoury and Weltman [10] and Gubser and Khoury [11] have shown that inclu- 
sion of non quadratic terms in the potential for a scalar field can greatly alter the experimental 
constraints on its coupling and mass. In particular, the range and effective strength of such a 
"chameleon" force can drastically depend on the environment. Such effects have been previ- 
ously explored in theories of time varying alpha [12-14]. Several subsequent works explored 
such terms in quintessence dark energy models [15-19]. In this paper, we further explore 
the experimental implications of the chameleon effect. We consider the constraints on the 
couplings to ordinary matter of a nearly massless chameleon field with WEP preserving cou- 
plings. Our results are dramatic: we find that the inclusion of an extremely tiny quartic 
coupling constant for the scalar field, as small as 10 -53 , weakens the constraints on the al- 
lowed coupling of the field to ordinary matter. We find allowed parameter regions for new 
forces which have not been discovered previously. Similar conclusions would apply to EP 
constraints, although we do not do a numerical analysis of those constraints here. 

2. The Model 

We consider a real scalar field theory governed by the Lagrangian density 

-£ = i(<90) 2 + vm (2.1) 

We interpret V{4>) as the renormalized effective potential. We will consider the (time- 
independent) solutions of the equations of motion from this Lagrangian. Our analysis uses 
classical field theory, which is valid for a weakly coupled quantum theory. We also assume 
that (j) has been suitably shifted by a constant so <\> = is a global minimum of V, and that 
V is analytic in the vicinity of = 0, which is equivalent to assuming that there is no degree 
of freedom which becomes massless at = 0. 

If we expand V about = 0, the term linear in must vanish by minimization of V. 
Furthermore, we assume the minimum should be a global minimum, so that we do not have 
to consider whether the coupling to matter could trigger formation of a catastrophic bubble 
causing tunneling to the true vacuum. Therefore, the size of the coefficient of the cubic term 
is bounded in terms of the coefficients of the quartic and higher order terms. For arbitrarily 
weak fields produced by coupling to small source terms it suffices to consider only the 
quadratic term. Macroscopic sources produce larger field strengths for which the cubic and 
higher order terms also become important. We will stop our expansion at quartic order. 
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Because we have assumed = is a global minimum of the potential, for a vacuum potential 
of the form 

V = Im 2 ^ 2 + ^ + -/ (2.2) 

the coefficient of the cubic term satisfies 

\p\ < 3^m 2 e/2 . (2.3) 

The contribution of the cubic term is then larger than that of the quadratic term only for 
field strengths large enough so that the quartic term is of at least comparable size. Thus 
when computing the effects of the force due to (ft, provided the quadratic and quartic terms 
are included, neglecting the cubic terms will lead to errors that are of most of order one. 

We assume <p couples to the fields of the standard model. Inside most ordinary forms of 
macroscopic matter the expectation value of any scalar operator involving the standard model 
gluon, quark or electron fields is approximately proportional to the mass density. To estimate 
at leading order the modification of the constraints on a scalar field due to the inclusion of 
nonquadratic terms in the potential it therefore suffices to consider the effective potential 

V{<t>) = \m 2 4> 2 + - Ptp, (2.4) 
where p is the mass density of matter. The equation of motion for </> is then 

d 2 cj) = m 2 (l) + €(p 3 - (3p = ^. (2.5) 

dcp 

In the limit e^Owe recover the Yukawa theory which is generally the form assumed for new 
forces. However in general for a scalar field there is no principle requiring that e = 0. In fact, 
a nonzero e is required to renormalize the theory. 

It follows from the Lagrangian density (2.1) that a test mass m! would experience a fifth 
force 

F$ = -(3m'V(f>. (2.6) 

Experimental tests of the gravitational inverse square law generally put constraints on 
the parameters a and A of an additional Yukawa term 

F = G N m 1 m 2 V . (2.7) 

r 

These parameters are related to those of eq. (2.5) by 



m= 1/A,= y/4irG N a (2.8) 

in the case e = 0. 

A Yukawa force falls off exponentially fast at distances larger than A = 1/m. The e(f> 3 in 
the equation of motion will cause (f> to fall even faster, at least as fast as a Yukawa field with a 
mass of m e s = \f m 2 + e0 2 . Following ref. [10], we refer to as a chameleon field because the 
associated force has an effective range and strength that varies according to the distribution of 
matter. In the next sections we show that the constraints on m and (3 are generally weakened 
by nonzero e, and present both analytic estimates and numerical computations of the effect. 
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3. Analytical approximations 



3.1 Inside a Large Object: 'thin shell' approximation 

The effective range of the force between test bodies is set by the inverse square root of the 
curvature of the potential. The chameleon potential is higher order than quadratic, and so 
the curvature of the potential about a nonzero field value may be much larger than it is about 
zero. When the range becomes shorter than the size of the object sourcing the field, then 
only a "thin shell" of the object [10] of thickness equal to the effective range of the force acts 
as a source for the field. 

For sufficiently large e, we may use a thin-shell approximation to find the value of <f> 
inside a large object whose density varies slowly relative to the effective range of the force. 
As a result, 4> inside the object varies slowly and remains always near the value 4> ma x that 
maximizes 

-V(4>) = ~m 2 4> 2 - + 04>P- (3.1) 
One condition required for this this approximation to be valid is: 

m eS = ^/V'% max = V™ 2 + 3eC a . > j (3-2) 
where i is the scale over which the source density varies 

P 

Note i < R, the radius of the source. The left hand side of (3.2) is an effective mass for the 
(p field, which depends on the density p of matter through the maximization of eq. (3.1). 

This approximation to minimizing the total energy neglects the energy of the gradient 
terms in the field. Since the field must go to zero far from the source, the gradient energy out- 
side the object may not be negligible. A second condition for the validity of the approximation 
is 

<t>max < 4>y{R)- (3.3) 

with 4>y{R) ~ — the Yukawa potential at the surface, and M is the total mass of the 
source. This condition ensures that, as e — > 0, we recover the Yukawa case. 

Maximization of (3.1) determines the overall strength of the field produced by a source 
with a thin shell. When 

£>> U) ■ 

the mass term in V{4>) is negligible, and (f) ma x is given by 

f (3p \ 1//3 

(pmax ~ ( — j • (3.4) 

Substitution of eq. (3.4) into the thin shell conditions (3.2) and (3.3) shows that the 
earth has a thin shell for ae^ 10~ 65 . 
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3.2 Near surface behavior of a thin-shell object 

For a spherical source, static solutions to (2.5) obey the ODE 

r a<p 

Exterior to a spherical object we can get a good understanding of the solutions to (3.5) by 
considering the relative importance of the nonlinear term. Since this term increases m e s, we 
expect (ft to decay fast when (3.5) is nonlinear, and to approximate Yukawa behavior when 
the nonlinear term is unimportant. Let us define 

71 ~ |m 2 0| + |0'/r|' { ' 

as a measure of the importance of nonlinearity in (3.5). 

When the equation is sufficiently linear, the solution has Yukawa behavior, (p' Hn « — (m + 
l/r)(pii n . For r <C 1/m, we can self-consistently approximate 

9lin ~ (3.7) 

setting r]i in < 1/5. If we attempted to apply this approximation right at the surface of a thin 
shell object, it would break down: 

, rD m cS R m efi e 1 

W/ ~ PmaxVtR r > > (3-8) 

by (3.2). So the conditions for a thin shell guarantee the nonlinear term is important near 
the surface. 

When the equation is sufficiently nonlinear, we can neglect the terms in the denominator 
of eq. (3.6). The LHS of (3.8) grows like (ae) 1 / 6 , so this is justified near earth's surface for 
most of our parameter space. Then our differential equation is 

<f>" ~ e</> 3 , 

which would be exact for planar geometry and m = 0. This equation has the monotonically 
decaying or growing solutions 



Matching to (p m ax at r = R gives 



± = . (3.9) 

^l(r-R± Jl^ox) 



So (f) decays extemely fast, behaving like it has a pole y/6 of a thin shell distance (i.e., m~^) 
inside the surface. Eventually, rj falls to a value where the linear terms become important. 
Neglecting m and using (3.9) in (3.6), one finds this happens at 

R - y/6/mrff 
riin ~ 1 - 2/r/ 2 ' 
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For a large thin shelled object, ru n is of order R. Even though (t> ma x increases with ae, the 
shell thickness decreases, and ru n is insensitive to ae. For instance, as long as ae> 10~ 52 , 
m c g of the earth is large enough that ru n is independent of ae to 1%. 

Clearly, <fi must fall fast as the shell gets thin. One may wonder if this makes the 
chameleon force easy to detect near the earth's surface. To answer this, we also use eq. (3.9) 
to derive the field gradient near the surface: 

This shows that even though 4>' (Re) / 4>(Re) grows with e, (p m ax falls quickly enough so the 
surface force always decreases with e. 

The picture we have derived here is that <p decays very fast until the linear terms in the 
ODE become important, then takes on Yukawa behavior. The net effect of e / weakens the 
effective coupling [3. All these results are confirmed by numerical simulations, although the 
location of linearization is not totally clear-cut. 1 

3.3 Effective coupling approximation 

The results of the previous section show that for a sufficiently massive thin shelled object 
coupled to a chameleon field, at a distance r which is large compared to the radius of the 
object but smaller than 1/m, the field strength is given by 

4>lin ~ "4=, (3.10) 

where c is a number of order 1. We can show that this result is always a good approximation 
for a spherical source whenever 

e> * ■ (3.11) 

~ aM 2 G N 

and m is sufficiently small, and for any geometry at a distance which is much larger than 
the size of the source and shorter than 1/m. Outside the source, for e = the approximate 
solution to the field equations is 

-mr 

(f> Y (r) = PVGnM . (3.12) 

r 

For nonzero e, (j) w iU De smaller, declining at least as quickly as a Yukawa field with mass 
of order m e ff. Thus if is larger than 4>u n , (j) w ih fell exponentially faster than 1/r until it 
decreases to a value of order <j>i in . 

Note that the differential equation (3.5) is covariant under the transformation 

<f>^kcj),e^ e/k 2 ,(3^ k(3 , 



x It is difficult to describe analytically the region between r) ~ 2, where linear terms become important, and 
r) ~ 1/5, where they become the only important terms. Numerical solution of the equations for m = shows 
that <j> scales like 4>u n oc e _1//2 instead of (pmax oc e -1 ^ 3 , for r very near the surface. But the r for which <j> 
scales with e~ 1//2 is smaller than the value of r for which 4> decreases with r as 1/r. 
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where k is an arbitrary scaling factor, which is consistent with the result that 4> scales as 
e-V2. 

Comparing (pn n with the Yukawa field that would be produced for e = we see that in 
comparing the chameleon field strength with a Yukawa potential, at a distance which is large 
compared with the size of a source one should use the effective coupling 



A. ~ (3.13) 
This approximation is good whenever the equations 

m<^l/R (3.14) 

and (3.11) are satisfied. 

Note that this effective coupling is independent of (3, and are effects of a chameleon force 
are weaker than gravity for any object which is much more massive than the Planck mass 
divided by y/e. Note also for any given object that the effective coupling depends only on e 
and the mass of the object. 

Thus if we consider any given experimental test of the inverse square law involving two 
sources of mass M\ and M 2 whose separation is larger than their size, which bounds a certain 
a and m for a Yukawa field, we conclude that for a chameleon field with the same m we have 
a lower bound on e which is 

1 

6 > M l M 2 G N a • (3 ' 15) 

Provided the inequality 3.15 is satisfied the chameleon is not constrained by the experiment 
in question. Similar conclusions would apply to equivalence principle tests, provided a is 
replaced by the product of a and the appropriate equivalence principle violating parameter. 



4. Experimental Constraints 

This section derives the major constraints on chameleon forces from different experimental 
searches for fifth forces. Again, these experiments generally put constraints on forces of the 
form (2.7) for different length scales A = 1/m. Our calculations of the modified constraints 
for e 7^ are summarized in the figure. 

Typically, a Yukawa force is constrained only for 1/m of order the characteristic lengths 
in the experiment. If m is too large, a Yukawa force decays before it is detectable. If m 
is too small, the Yukawa force appears to obey the inverse-square law. Here we consider 
how e / modifies this behavior (always weakening the force, as discussed above). But we 
also consider how a modification to Gn might be detected by affecting different experimental 
measurements of Gjy differently. 
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4.1 Lunar Laser Ranging 

The Lunar Laser Ranging experiment [20] puts the tightest constraints of any experiment on 
a Yukawa force, and the fact that it is loosened for extremely small values of e is our most 
prominent result. For this reason, we chose m « 1/ REarth-Moon, the most constrained value 
of m when e = 0. 

The experiment carefully measures irregularities in the moon's orbit by reflecting a laser 
beam from a reflector on the moon. This produces a constraint on Yukawa forces of a ^ 10~ n 
for tyi ~ 1 / 1 ^Earth—Moon [^]* Let us apply the effective coupling approximation derived in 
section 3.3. We need to determine a e ff, the strength of the chameleon force on the earth- 
moon system relative to the strength of a Yukawa (e = 0) force with a = 1. 

We start by assuming both the earth (e^ lCT 65 /a) and moon (e^ lCT 61 /a) have thin 
shells. The chameleon force is given by 

Fc = V REarth _ Moon J d*x l -{Vcpf + V(cf>), (4.1) 

where (ft is a static solution to eq. (2.5) with both the earth and moon as sources. As we saw 
above, a thin-shell source produces a field that falls to 4>iin in a distance roughly ~ Rsource/ (1 — 
2/rjf in ) independent of e. Farther from either source, the behavior is essentially linear, so (f) ~ 
4>Earth + 4>Moon with (j) Earth and (pMoon single-source solutions, since the nonlinear regions do 
not overlap. Since most of the overlap occurs in the linear region, one can by definition neglect 
the €(p s term in V((p), and the force too behaves linearly. So it is separately proportional to 
/? e ff,i, the single source apparent strength of each field: 

Arfr,E/?eff,M 
acS = 4nG N ■ 

After reaching 4>u n , the single-source field has the Yukawa form, so 



a n 4>lin ^Vlin 

where we have used e~ mr ' in ~ 1. Thus 

UcS ~ eG N M E M M (42) 
is the effective strength of the chameleon force (the e —> limit clearly cannot be taken once 
a thin shell approximation is used). Our discussion above anticipates that a e g is insensitive 
to a, since the large r values of the field are insensitive to the surface value, which by eq. 
(3.4) is (ftmax oc a 1 / 6 . 2 Note, though, that the thin-shell conditions guarantee a e g < a. Eq. 
(3.3) in spherical geometry gives 

3 ( 4 ! 
e > 7T 



PM 



2 A 10-fold rescaling of the surface value, equivalent to a rescaling of (3 by 10 3 or a by 10 6 , barely affects 
the later field. 
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hence 

r?L min(Mi,M 2 ) 

< «— k^FT < a / 100 - 

7r max(Mi,M 2 ) 

Then the LLR constraint a c g < 1CT 11 gives: 

47r1D 11 r7 2 

e > ^ „ 5 . 10-53 ae > 10-6I. (4.3) 

G N M E M M ' ~ V ; 

when both sources have thin shells. 

For ae intermediate between the earth's and moon's thin shell values, we treat the moon 
as a point test mass. Then 

PPef£,E I a n 
«eff = . ' = W . ^ Asff,E < 10 , 



47tGat V 47tGat 
so 

in 44 

« < ^77—^ ~ e(2 • 10 43 ), 10" 65 < ae < 1(T 61 . (4.4) 

As ae lowers to the value for which the earth acquires a thick shell, this constraint neatly 
approaches 10 -11 , the e = constraint. 

4.2 Helioseismology 

One way to phrase the gravitational inverse-square law is to claim that the gravitational 
constant Gn is indeed a constant. So comparisons of measurements of Gn with objects of 
greatly different size can constrain our model if the chameleon causes an apparent modification 
of Gjy. In section 4.3, we will consider the constraints from comparing terrestrial Cavendish 
experiments, which involve source masses from ~ 10 kg to ~ 10 4 kg. But first, we compare 
this whole group of experiments with Gn determined by data from the sun, which has a mass 
~ 10 30 kg. 

The standard solar model balances pressure gradients against gravitation, fixing the grav- 
itational force between volume elements when the pressure is known. Solar neutrino data 
constrains the temperature, and thus the gas pressure, to about 1% [21]. This constrains Gjy, 
and not just the product MGn, to 1%, in agreement with lab-based Cavendish experiments. 

^,From eqs. (3.3) and (3.2), we find that the sun has a thin shell for ae^ 1CT 75 . A thin 
shell means that (f> is constant inside the source (or varies quasi-statically for slowly changing 
density), so the chameleon force is then negligible in the sun's interior. (Strictly speaking, for 
thin shell we should compare the effective mass to the length scale £ for density variations, but 
we approximate by comparing to the solar radius). The largest-scale Cavendish experiment 
has a thin shell for ae ^ 10 -23 . 

In this and the following sections, we derive strict constraints and simplify our analyses 
by assuming the chameleon force completely negligible in each experiment when the source 
has a thin shell, and completely Yukawa when thin-shell fails. This simple picture is valid 
except for ae near the endpoints of the thin shell regime, when r\ ~ 1 and the nonlinear and 
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linear terms in the ODE are of the same order. Improving the approximation would weaken 
the constraints in this regime, smoothing out their sharp edges. 

In this regime the Cavendish experiments see a Gn modification of relative strength a, 
while the sun sees only gravity. Thus, we derive the constraint: 

For 1(T 75 <ae< 1(T 23 , a < 1(T 2 . (4.5) 

4.3 Lab Cavendish experiments 

The modern Cavendish experiments [22-32] involve a range of techniques, geometries, and 
source masses from ~ 10 kg to ~ 10 4 kg. There is some disagreement among results, and 
it is amusing to note that the chief outlier [23] differs from the accepted value by the right 
sign to be explained by a chameleon. However, to derive the most stringent constraints, we 
discount this outlier and assume the other experiments accurate to their stated uncertainties. 
The experiment with the smallest source masses [31] used 4 stainless steel spheres of mass 
8 kg, giving a thin shell for ae> 10~ 17 . The largest source mass was [32] a pair of mercury 
tanks 7550 kg each, with a thin shell for ae^ 10~ 23 . 

Then a must be such to make these two experiments agree within their fractional uncer- 
tainties of 10~ 5 : 

For 10~ 23 <ae< 10 -17 , a < 10~ 5 . (4.6) 

4.4 Ocean and Lake Experiments 

Zumberge et al [33] measured g at varying depths in the ocean as a test of ISL and a de- 
termination of Gat. Their value agrees with lab Cavendish experiments to about 3 • 10~ 3 . 
Again, we expect no chameleon force in the ocean when the earth has a thin shell because 
4> is quasi-static in the interior. This gives a constraint when the earth has a thin shell but 
laboratory experiments do not: 

For 10- 65 <ae<l(T 23 , a < 3 • 10 -3 . (4.7) 

The lake experiments [4] use a lake as a source and typically find a value of Gm that 
agrees with lab values to about 10~ 3 . This constraint is essentially subsumed by the ocean 
experiment. 



4.5 Tower Gravity Experiments 

Several experiments [4] measured constraints over scales of a few hundred meters by detecting 
deviations from an inverse-square falloff of g as one ascends a tower. This is the case, near 
the surface of the earth, where the chameleon is strongest and falling fastest, among all the 
experiments. 

They observed limits on deviations of ~ 10~ 7 <7 = 7 • 10~ 17 /Re, setting h = c = 1. For 
the earth with a thin shell, we use our approximation eq. (3.9) to compare the acceleration 
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field a distance ~ 10 2 m ~ 10 4 Re = X^e above the surface to the expected ISL scaling: 



(3 



(1 + X? ' 

7- 1(T 17 
Re 



1 



< 



(X + e - 1/6 V2(3^ i ) 1 /3)2 (1 + X) 2 26-V3(^_)2/3 

(4.8) 



Substituting in numbers: 



(10~ 4 + (ae)' 1 /6(2.2 • 10" 11 )) 2 (ea)" 1 /3(5 . 10" 22 ) 



< 5 • 10 24 



(4.9) 



For ae <C 10 40 , the shell becomes much thicker than X&E, and the constraint approaches 



a 



1 



1 



(ae)V6 2.2 • 1Q- 11 



<10 7 , 



which becomes irrelevant for ae ~ 10 -64 , just as thin shell breaks down. For ae S> 10~ 40 , 
the first fraction in eq. (4.9) can be neglected, and a^S 10 4 e 1 / 5 . For 10~ 21 , the overall 
chameleon strength is small enough that the constraint weakens to 1. 

4.6 Lab Inverse- Square Law: Spero et al 

Some of the tightest lab constraints on Yukawa forces in the centimeter length scale come 
from Spero et al. [34,35]. This null experiment involved placing a 20 g test mass about 1.2 cm 
from the central axis of an Fe cylinder of inner radius ri n = 3 cm and outer radius r ou t = 4 
cm. The experiment used a cancellation mass to cancel edge effects, so the force should vanish 
inside the cylinder for inverse-square law. Spero et al. observed a torque corresponding to an 
acceleration field less than about 3 • 10 -12 m/s 2 . 

As usual, our field becomes linear for a thick shell, so constraints are only different from 
those on a Yukawa field when e large enough to give a thin shell, but small enough that the 
force is still detectable at the test mass position. The thin shell criterion (3.2) gives 

VsiPp) 1 ^ 6 > — 1— , 

Tout fin 

so the Yukawa constraints are modified when ae^ 1.5 • 10~ 14 . Notice that a large value of e 
is needed to reduce the shell distance to a centimeter. 

The cylindrical geometry changes 2<p ! /r in eq. (3.5) to (f>'/r, but we neglect this term 
anyway in our usual approximation (3.9) for very nonlinear fields. 3 Then the acceleration 
field at r = 1.2 cm is 



/ty'(2r in /5) 



e/2/3 



<a o6s «l.l-10- 29 /r 



(V^r m (l -§) + (jfc)V3)2 



The validity of this approximation should extend 



a few cm, or most of the way inside the 



l-„2/2 

cylinder. Numerical results confirm the approximation, at least to order of magnitude, about 80% of the way 
from the wall to the axis. 
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This gives the constraint 



a 



■<e 1 / 3 (2.5-10 6 ). (4.10) 



( e l/6 + a -l/6 6 . 6 . 10-3)4 

This constraint equation can be solved numerically for a given e satisfying thin shell. The 
constraint would weaken even further for e ^ 3 • 10 -10 , as the test mass acquired a thin shell. 

Note that eq. (4.10) depends sensitively on the position of the test mass. Changing the 
position of the test mass would change the factor 2/5 in the first equation above, which would 
introduce a factor changing the relative importance of the two terms in the denominator of 
eq. (4.10). As the position of the test mass approached the inner wall, the constraint would 
strengthen to a^$e 1//5 (4 • 10~ 2 ). However, to our knowledge no such experiment has been 
performed. 

4.7 Lab Inverse- Square Law: Hoyle et al 

The University of Washington Eotwash group [36] put tight constraints on a Yukawa force 
in the 100 micron range. The technique was to measure the torque on a torsion balance 
pendulum produced between a pair of discs with holes bored in them. In this setup, the 
"missing mass" of the holes is the attractor for any field obeying a linear differential equation. 
But it is difficult to model such a problem for a nonlinear field. 

The basic premise of this null experiment is that the force between large parallel plates 
is independent of separation for Newtonian gravity. The setup made use of three discs, the 
lower two of which (7.8 mm and 1.8 mm thick Cu) were at a constant separation and had 
the holes offset azimuthally by 18 degrees. The upper disc (2.0 mm thick Al) was at a 
variable separation of 0.2 mm to 11 mm and acted as a torsion pendulum. Because the two 
lower discs were azimuthally rotated, they largely canceled each other's Newtonian torques 
on the pendulum. But for a short-range force, only the upper disc would exert torque on the 
pendulum. So we get to the core of the problem by calculating the separation dependence of 
the force between parallel plates for the chameleon field. 4 

The chameleon force becomes short range when linearity fails, which happens roughly 
when the source acquires a thin shell. The two lower discs were flush, so we can consider 
them as a single attractor with holes partway through. Then the experiment measured torque 
caused by the differential force from the parts of the attractor closer and farther from the 
pendulum. For ae< 6 • 10~ 10 , neither pendulum has a thin shell, and for sufficiently small 
m we expect no ISL deviation in the experiment. For e very large, we expect the force 
to be negligible. For intermediate values of e, we derive the constraints by computing the 
ratio of the forces from the two attractor discs and comparing to the Newtonian value of 
approximately 1 within the experimental sensitivity of about 1%. 

For this planar one dimensional problem, our ODE becomes 



4 Of course, the cancellation of torques was not exact at all separations, but there is more cancellation for 
linear than for nonlinear fields. 
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with the origin chosen halfway between the plates. It follows from the differential equation 
that the quantity 

C = <f>' 2 - 2V(<f>) (4.11) 

is independent of position, in analogy to conservation of energy in a 1-d mechanics problem. 
The potential energy density in the field is given by 

u = cj)' 2 /2 + V{(f>) = C/2 + 2V(<f>). 

Then the 1-dimensional chameleon "force" between plates is the derivative of total potential 
energy with respect to separation: 

dU d f Ax ' 2 , 1, /A lns , . / \ \ f Ax/2 du , 

F c = — — — = — — — / udx = — (u (Ax 2) + u(-Ax/2)) - / tt^— cfe. 
d(Ax) d(Ax)J„ Ax/2 2 y y 1 ' y 1,1 J„ Ax/2 dAx 

So far, this is exact, but integrating (4.11) and inverting for is & n impractical way to 
find F c . If we apply thin shell, 0(-Ax/2) = 4> max ,Cu, 4>{-Ax/2) = 4> max ,Ah then u{±Ax/2) 
depends on separation only through C/2. We further assume that, for our parameters, does 
not change much from (j) max - Then it is valid to approximate by a series solution to second 
order: 

A ~ A , A i / 2 i i ( §max,Cu ~ <Prnax,Al \ ( £ ,3 m 2 \ 2 

pa 0o + 4>ix + ^22; =0o+l ^ — J x + I + "Y^o J z 

and match to the thin shell solution: 

with ro = 0.1 mm. The last equality can be solved numerically for 0o as a function of e and 
a. Then 

C = (0i + 20 2 a;) 2 - 2^(00 + fax + 2 x 2 ) 
gives us the constraints. We compare to the Newtonian 1-force 

itGn , . 

= 2~ (P2*2 - Pl*l) , 

with ti the thickness of plate i, by comparing 

F c (Ax) + Fjv 



F C (Ax + t upper ) + F N ' 



with t upper = 1-8 mm, to 1. 

We expect from (4.11) the ODE will become linear when V{4>) ~ V{<j) m ax) *C C. This 
occurs for at <C 4 • 10~ 4 for separations of 0.2 mm and ae <C 3 • 10~ 14 for separations of 10 
mm. This latter inequality gives us some idea of the constraints since a linearized field is 
likely to behave like the Newtonian case. 
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Indeed, numerical simulations show that 7 • 10~ n /e for 10~ 8 . For larger e, the 
chameleon force weakens and the constraints become less relevant. For e = 10~ 6 , we have 
a ?S 0.1, and for e = 10~ 4 , a ?S 1. Beyond this, the parabolic approximation breaks down, but 
it is safe to assume that a ^ 1 is allowed. 

Hoyle et al also did a calibration using metal spheres of radii 4 mm and 2.5 cm, separated 
by 14 cm. For a ^$ 10~ 15 /e, all the calibration spheres would have thick shells, so the measured 
force is 1+a times the gravitational force for the calibration spheres as well as for the pendula. 
There is therefore no constraint on a in this regime. For 10 -15 ^ ae^S 6 • 10 -10 , the larger 
spheres but not the pendula would have a thin shell. The chameleon then contributes more 
strongly to the pendulum force than to the calibration force. Since the precision is about 1% 
and no deviation is observed then a must be less than 10~ 2 in this regime. 

4.8 Lab Inverse- Square Law: Hoskins et al 

Hoskins et al [35] put constraints on the ISL at the 10 cm length scale. They compared the 
force on a torsion pendulum from a 7 kg "far" mass at about 105 cm separation to that from 
a "near" mass at 5 cm and found a^z 10~ 3 . However, their experiment does not appear to 
constrain a chameleon force. Since both masses were copper, they would have radii about 5 cm 
and 1 cm respectively. But because each mass was located at least ~ 5 radii from the torsion 
pendulum, one would expect that the field from each mass at the torsion pendulum would 
always be linearized, either because each source had a thick shell or because the field from a 
thin shell source had already decayed past 4>u n - Since we use small m, a linear <\> field would 
look like Newtonian gravity in this experiment. We have already considered modification to 
Gn in Section 4.3. 

4.9 Planetary 

Planetary constraints [4] are negligible at distance scales large compared to 1/m even for 
e = 0, and as we have noted, the force always weakens with e. If we considered the 1/m ssAU 
case, the discussion would parallel that for LLR. 

4.10 Free Fall 

It is interesting to consider that, in the regime where two objects both have thin shells, the 
chameleon force behaves like a Yukawa interaction with a c g- independent of a (see Sections 
3.3 and 4.1). This is true when the separation is larger than ru n ~ R SO urce/ (1 — 2/r]f in ). If we 
consider the regime where satellites have a thin shell (ae^ 10~ 18 for a 10 kg metal sphere), 
then satellites at heights greater than ~ 5Re would experience a fifth force with a e g- given 
by eq. (4.2) dependent on mass, affecting their periodicity. 

From the discussion in Section 4.1, a c g would be down from a by a factor M sa t/ME ~ 
10~ 20 even for the minimal value of e which gives a thin shell, and would decrease with e. So 
such a fifth force would be negligible. 
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Figure 1: Plot of allowed region for the coupling a, as a function of the nonlinear term e, with the 
mass of the field set to the inverse earth-moon distance. The constraints do not change appreciably 
for massless or for somewhat more massive fields, except in the region e < 10~ 50 . This entire region is 
forbidden by LLR when e=0. 

5. Conclusion 

Currently, there is no evidence for deviations from Newtonian gravity, and this concordance 
is generally interpreted in terms of constraints on new forces. Usually, the equations for any 
new force field are assumed to be linear. In this work, we have examined how such constraints 
would be reinterpreted in the presence of a nonlinear, "chameleon" force. We note that some 
constraints on the strength of such a force become much weaker, even for extremely small 
nonlinear terms, and there is room to hide a new long range force whose effects could be as 
large as 1% of gravity between small objects. 

Although currently there is no compelling outlier, it is worth emphasizing that for 
chameleon forces, the possibility arises for different experiments to produce inconsistent re- 
sults, when interpreted in terms of the Yukawa framework. Therefore outlying results on 
searches for new forces should not automatically be dismissed without further investigtion. 
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